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Abstract 

We analyze the long time behavior of initial value problems that model a process where particles of type A and B 
diffuse in some substratum and react according to nA + nB — * C. The case n = 1 has been studied before; it presents 
nontrivial behavior on the reactive scale only. In this paper we discuss in detail the cases n > 3, and prove that they show 
nontrivial behavior on the reactive and the diffusive length scale. 
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1 Introduction and main results 

We consider, for arbitrary but fixed n G N, ti > 1, the reaction-diffusion problem 



bt = bxx ' 

for a; G R, f > T > 0, with initial conditions a(x,T) = ao(a;), b{x,T) — bo(x), satisfying 

Um ao(a;) — 1 , 
= 1 , 

= Km bQ{x) — . 



(1) 
(2) 



lim 6o(x) 

X — * + oo 



(3) 



lim ao{x) 

X — >+oo 



The choice of the initial time t ~ t, and a class of initial conditions oq, 6o will be described later on, but for the purpose 
of this introduction it is useful to have in mind the "natural" case: r = 0, ao(x) = 1 for a; < 0, ao{x) — for x > 0, and 
bo{x) = 1 for a; > 0, bo{x) = for x < 0. 

This initial value problem models the time evolution of a chemical system of two (initially separated) substances A and 
B, that diffuse in some substratum and react according to nA + nB — > C, with a substance C that is supposed not to 
participate in the reaction anymore. The model is a mean-field description of such a situation where the functions a and b 
represent the densities of the substances A and B. For more details see [||. 

Equations (0) and ^ are best studied in terms of the sum 



V = a + b 



(4) 



and the difference 



u — a — b , 



(5) 



which satisfy the equations 



ut = u^x , (6) 
vt = - {v^ - v?T , (7) 



with initial conditions vq and uq (at time t = t) that satisfy 

Km uq{x) — 1 



and 



Km uq{x) — —1 , (8) 



Km v^(x) = 1 . (9) 

X — *±oo 



For initial conditions ao, bo with 

ao(x) = hn^-x) , (10) 

the functions vq and uq are even and odd, respectively, and the equations (^, (^ preserve this symmetry. Furthermore, for 
the special initial condition 



u{x,t) = -fIj^{x/^/T) , (11) 



with ijii defined by the equation 



/^i(y) = crf(|)^--= / e-'^da, (12) 



equation (^ has the explicit solution 

u(x,t) — —fii{x/Vi) . (13) 

We note that the initial condition (^ij) for u (at time i = r) is simply the solution of equation (||) with the "natural" initial 
condition, u{x,0) = 1 for a; < 0, u{x,0) = —1 for a; > 0, evaluated at < = r. To keep this paper as simple as possible we 
now restrict the discussion to this case, i.e., we consider from now on equation (^ with initial conditions satisfying (^, and 
u given by ([T^). We note, however, that more general (asymmetric) initial conditions for u could be treated as well. This 
would lead to corrections to u of the order 0{l/t), and such corrections do not change in any way the discussion of the 
equation for v that follows. 

The reaction-diffusion problems considered here develop, in addition to the built-in diffusive length scale 0{^/t), an 
additional shorter length scale, on which the reaction takes place. The function F, 

F^i(4a6)"^i(«2-^2)", (14) 

is called the reaction term or reaction front, and we are interested in describing the asymptotic behavior of the function 
F for large times. The knowledge of this behavior is useful, since it appears to be universal, in the sense that it is largely 
independent of the choice of the initial conditions and of the details of the model under consideration. As mentioned above, 
if Vq is an even function, then v and as a consequence F are even functions of x. We will see that the critical point of F at 
X = is a maximum, and that F decays (rapidly) for large x. 

Before proceeding any further we note that the factor of 4"~^/^ in (|l]), (||) and (|l^ ) is just a normalization, and has been 
chosen for convenience to make the equation (0) for v look simple. In fact, any system of the form 

at = Dattxx - ka^ahY , 
ht = 0^3,^,^ - U(ahY , 

with positive Da^ Dh, ka, and k^, and with initial conditions such that 

lim a(x, 0) = aoc > , 

X — *■ — oo 

lim b{x, 0) = feoo > , 

X — >oo 

lim a{x, 0) = lim b{x, 0) = , 

X — > + C30 X — ^ — OO 



can be reduced, by scaling space and time and the amplitudes, to the problem 



at = axx - ^(4a6)" , 

with I? > 0, and with initial conditions such that 

Hm a{x, 0) = 1 , 

X — * — oo 

hm 6(a;,0) = /3 > , 

X — >oo 

hm a(x, 0) = hm 0) = . 

X — * + oo X — ^ — oo 

In this paper we have hmited the discussion to the case (3—1 and D — 1. The case j3 ^ 1 leads to a moving reaction 
front. A change of coordinates to a moving frame complicates the analysis, but the problem could still be treated with the 
methods presented here. Choosing D = 1 makes the equations mathematically simpler. As a consequence, as we have seen, 
the two equations for a and b can be reduced to just one equation for the sum v — a + b, since the equation for the difference 
u — a — b can be solved explicitly. Even though we do not expect the asymptotic behavior of the solution to change in any 
relevant way ii D ^ 1, the strategy of proof would have to be changed considerably, since the equations can not be decoupled 
anymore in that case. 

Before we state our results, we briefly discuss the expected dependence of the results on the parameter n. 
The case n — 1 has been studied in detail in where it is proved that in this case the reaction term ( p^ satisfies, for 
all z g R, 



lim t^'<F{ez,t) = p{\z\) , 



where a = 1/6, and 7 = 1/3, and where p: R+ is a smooth function that decays like exp(— const. z'^/^) for large values 

of z. It follows furthermore from the results in that the function F is very small on the diffusive scale in the sense that 
for n = 1, y 7^ 0, and all p > 0, 

lim iPF(Viy, t) = . (15) 

t — >oo 

The smallness of F on the diffusive scale is easily understood by realizing that, for n — \ and for positive values of x on the 
diffusive scale, i.e., for x/\/t >> 1, equation (|l|) essentially reduces to 

at = axx ~ Aa , (16) 

with A > 0. Therefore, the function a decays exponentially fast to zero on this scale, and similarly for b for negative values 
of X. 

For n > 1, however, equation (|l|) reduces, for x/^/t >> 1, essentially to 

at = axx ~ Aa" , (17) 

with A > 0. The solution of ( |l7| ) has an asymptotic behavior that is radically different from the solution of (|l^). In 
particular, for n = 2, the solution may even blow up in finite time if a is not a positive function. Note that, for n odd, the 
nonlinear term in (|l^ is always a "friction term", independent of the sign of a, and the case of n odd will therefore turn 
out to be easier to treat than the case of n even. It is well known |^ that for n > 3 and small bounded integrable initial 
conditions, the nonlinearity in ( p^ becomes irrelevant for large times in the sense that the solution converges to a multiple of 
exp(— a;^/4t)/v^, which solves the linear equation at = axx- We would therefore expect that, for n > 3, the function F is of 
the order 0{t~^/'^) on the diffusive scale. This turns out to be wrong. As we will prove below, F is of the order 0(t-"/("-i)) 
for n > 3, because F converges on this scale pointwise to a function that is not integrable at the origin. This corresponds 
to a solution of (^7|) for which the nonlinear term is a marginal perturbation, i.e., a solution with an amplitude of the order 
see that one can take advantage of this fact, and a diffusive stability bound will be good enough to 
prove convergence of F to its limit, but as a consequence, our results will be limited to the case n > 3. The cases n — 2 and 
n — Z are special and will not be discussed any further. 
The following theorem is our main result. 

Theorem 1 For arbitrary but fixed n G N, n > 4, there exist t > 0, functions pi /Xj, ifi, (^21 and a class of initial conditions 
(specified at t — t), such that (j^ has a unique solution v that satisfies for all t > t the bound 

\v{x,t)~Vooix,t)\<^^^ , (18) 



whe 



v^x, t) = Mi(M) + t~^p,C-^) + t->i(^) + , (19) 



Remark 2 This theorem is a local result, in the sense that the class of initial conditions will be a set of functions in a 
(small) neighborhood of the function Uoo.Oi ^^00,0(2^) = Vooix,T). In particular, our methods do not allow us to show that the 
solution with the "natural" initial condition vq = 1 at t = belongs to this set at t = t. We do expect, however, that this is 
the case, as has been proved for n — \ in W. 



Remark 3 We note that, if an initial condition Vq is such that vo(x) — \u{x/^/t)\ < for a certain x, then ao{x) < 0, 
if X > 0, or bo{x) < if x < 0. A priori, we do not need to consider such initial conditions, since in our model a and b 
represent particle densities, and the solutions a and b are positive if the initial conditions ao and bo are positive. As we will 
see, for n > A, it will not be necessary to impose that ao and bo be positive everywhere, and it will neither be necessary to 
impose that vo = ao + bo be an even function. 

As we will see, the functions (pi and are small on the diffusive scale, i.e., for x « Vty, y 0, and t large, 

VooiVty,t)^^^,{\y\)+t-'^l,{\y\) + 0{t-''') , (20) 

where e' ~ e if n > 5, and 27 < e' < e if n = 4, 5. Using the definition (^), (|^) for v and u, we therefore find that for > 
and t large, 

a{Viy,t) = ^t-'n2{y) + 0{t'^'') , 

and similarly for 6, for y < 0. In contrast to the case n — 1, where only exponentially few particles reach the diffusive scale, 
the amount of particles decays only slowly for n > 3. Our results imply that, for large times, the density of the remaining 
particles is given by the function /i2, i.e., it is independent of the initial conditions. 

As a corollary to Theorem |l| we get a precise description of the reaction front F on the reactive and the diffusive scale. 
This description will be given in Section 4, once we have defined the functions /i^, fi ^^'^ i'^ Section 3. In Section 
2 we explain our strategy for proving Theorem |^. This strategy is implemented in Section 5 and Section 6. The Appendix 
contains the proof of the existence of the functions fj,2, fi and ip2. 



2 Strategy of the Proof 

Consider functions v of the form 

v{x,t) Vca{x,t) +ip{x,t) , (21) 

with Woo as in Theorem]^, and ip{x, r) = iPq(x), for some r >> 1, with ipQ £ LiD Lao. Substituting ( pl| ) into (|^) leads to an 
equation for the function -0 of the form 

^ = ^" ~ Vi; - I ~ T{%l)) , (22) 

for certain functions V and /, and for T some nonlinear map. We will show that if Voo is defined correctly, r large enough and 
■00 small enough, then V can be chosen positive and T will be small, so that the solution of equation (^2|) will be bounded 
for large times by the corresponding solution of the inhomogeneous heat equation xjj = ip" — I. We will find that, with the 
right choice of Voa , 



dx 



I{Vtx,t) < const, ^ (-23) 



from which the bound (|T^) will follow. We note that 47 < ^ for 71 > 4 > |, so that contributions of initial conditions will 
become irrelevant for large times, i.e., the solution ip becomes what is called "slaved to the inhomogeneous term". 



3 Asymptotic Expansion 

In order to implement the strategy outlined in Section 2, we need a function Vao that approximates the solution v for large 
times sufficiently well, uniformly in x. Since we would like to control the time evolution of equation ( p2[ ) on Li H L^o, this 
function Vao needs to satisfy lim2,^±oo Woo(a;, <) 1 in order for v to satisfy the boundary conditions (^. Furthermore, 
the inhomogeneous term / in equation (|2^) contains second derivatives of Woo, and the function / can therefore only be in 
Li n Loo if Vcx) is at least twice differentiable. We now construct a function Uoo satisfying these requirements through a two 
length-scale asymptotic expansion. 

To simplify the notation later on we use the convention that, unless stated otherwise. 



_ x 



(24) 



and 



t° 



(25) 



and we will refer to y as the diffusive length scale and to z as the reactive length scale. 

The function Voo is given by the first and second order terms of a so called "matched asymptotic expansion" . The 
"matched" refers to the fact that such an expansion contains functions that can not be associated uniquely with one of the 
length scales and can therefore be used to "match" the behavior at large distances of the shorter scale with the behavior at 
short distances of the larger scale. Let 



1 ' 



(26) 



-5 — 



-^y2-<5 = t-^^z^~^ and 



and let 7, a, and e be as in Theorem Then, the functions y = f'^z, t~'^y' 
t~^'^y~'^^~^^ = t~^'^ z~'^^^^ are of this form and will naturally show up in the function Voo- As a consequence, the representation 
( p^ ) for Woo is not unique. If we choose (as we will) to compute the expansion for v^o in the order of decreasing amplitudes, 
i.e., if we first compute the term of order 0{t^"'), then the term of order 0{t^'^), and finally the term of order 0{t^^''), we 
get a representation of of the form 



v^{x,t) = ^,,i\y\) + t-''r^{\z\) + t-- (/i2(|y|) - A \y\-') + t^^-'MM) , 



(27) 



where 



viz) = Viiz) + 



-5 



(28) 



with A a certain constant to be determined later. 

We note that, by definition, Voo is a symmetric function, and it is therefore sufficient to consider positive values of x if 
we choose appropriate boundary conditions at x = to ensure regularity. 

Finally, since we will need to describe the asymptotic behavior of various functions near zero and infinity, we introduce 
the following notation. Let / be a continuous function from R_|_ to R, A; a positive integer and pi < P2 < ■ ■ ■ < Pk real 
numbers. Then, we say that 



near x 



0, if 



i=l 



and we say that 



near x = 00, if 



lim 

2:^0+ xPi" 



f{x)-Y^f,xP^ 



i=l 



lim x^^ 

x^oo 



f{x)-Y^f,x-P^ 



3.1 Equation for (i-^ 

To lowest order the function v is asymptotic to ^i{\y\), with p.^ as defined in We note that pi has near y = the 

expansion 



/^i(y) = Ky + Ksy^ + 
where k = and K3 ~ ^ j^^- Furthermore, lim^^oo /^i(y) = 1- 



(29) 



3.2 Equation for ip-y 

We make the ansatz v{x,t) — /ii(|y|) + 1 •'lydzD which we substitute into equation (0). We multiply the resulting equation 
with t''^"^", and take then the limit t ^ 00, keeping z fixed. This leads to the differential equation 



77" = (2KZ'n + jf)"' , 



(30) 



where z is now considered a variable in R+. Since lim^^oo l^iiv) — I7 the correct boundary condition for 77 at infinity is 

lim 7]{z) = , (31) 

Z — >-oo 

and at 2; = we impose 

V'{0) = , (32) 

which makes the function /i]^(|i/|) + t~'^ri{\z\) twice differentiable at x = 0, since fj-i and rj are twice differentiable at zero and 
dx (/^i(|y|) + '7(1-^1)) (0) = 0. A proof of the following proposition can be found in the appendix. 



Proposition 4 For n G N. n > 2, there exists a unique function rj: R that satisfies equation (3C) with the boundary 

conditions (pi) and (jS^). The function 77 is positive, and has near z = the expansion 



viz) = 'na- KZ + ri^z^ - rj^z^ + 
with positive coefficients rj^, 772 '^'^'^ "Hi- z large, rj is of the form 



for a certain constant Aqo, with S as in (ko), 



and 



S' = 



(y4^5i(5TT)TT - 1] /2 2<n<b 
25 + 1 71 > 6 . 



We note that 3 < S' < 26 + 1. The constants r/g, 772, 774 and Aoo are given in the appendix. The function ip^ is defined in 
(08) in terms of 77. 



3.3 Equation for ^2 

We make the ansatz v{x, t) — ^i(|y|) + 1"'''77(|2;|) + i~^(M2(|y|) ~ ^ \y\^^) which we substitute into equation (|^). We multiply 
the resulting equation with i^+'^j and take then the limit t — > 00, keeping y fixed. Since rj{z) = rj{t'^y) and limt^oo t'^~'^rj{z) — 
Xy~^ ~ 0, this leads to the differential equation for the function /i2j 

+ lyt^2 + = : (34) 

where y is now considered as a variable in R+ . At 7/ = we impose the boundary condition 

limu2(2/)2/' = A, (35) 

which removes the leading singularity of the function /^2(l2/l) ~ at ?/ = 0. As we will see, the sub-leading singularity is 

proportional to \y\^ , which is not a twice differentiable function a,t y — (except for 77 = 4 where S — 2). This singularity 
will be cancelled by imposing appropriate boundary conditions for the function 952- 

The correct choice of boundary conditions for ji2 at infinity is somewhat less obvious. In the appendix we show that the 
condition limy^oo A*2(y) = is not sufficient to single out a unique function jj,2- If M2 does converge to zero at infinity, then 
it is asymptotic to a solution of the equation 

ji" + ^yji' + eji^O . 

This linear equation is compatible with a (very slow) algebraic decay, ^2(2/) ~ or with a modified Gaussian decay, 

^2(2/) ~ exp(— ?/^/4)/7/^~^^, with the algebraic decay being the generic case. It will be essential in later sections that ji2 
decays rapidly at infinity, and we therefore impose the boundary condition 

lim ji2iy)y^' = . (36) 



A proof of the following proposition can be found in the appendix. 



Proposition 5 For all n > 4, there exists a unique positive function /ij : R+ R that satisfies equation ( p^ /) with the 
boundary conditions and (p^. For ?/ small, the function /ij is 0/ the form 

(y) = Ay-* + Aoy2-* + Aiy^-* + . . . , (37) 

with 

^ _ l\ -2nK:,5{6 + \)-K{5-2e) ^ „ 

^°~2^(n-l)5(^ + l) + 2(2<5-l) ^^^^ 
wzi/i Ai 7^ and with A as m Proposition 0. For ?/ Zarge, the function /ij decays rapidly in the sense that 

M2(y) = exp(-^)(^ + ...) , (39) 

for some constant C > 0. 
3.4 Equation for (/92 

We make the ansatz v{x,t) = fJ-iUvl) +t~'^r]{\z\) + i~'^(M2(|y|) ~ A |?/|"'') + t-'''''(y92(|z|) which we substitute into equation (0). 
We multiply the resulting equation with ^s^+aa^ g^^^j ^g^j^g -^iign the limit t ^ 00, keeping z fixed. This leads to the (linear) 
differential equation for ip2, 

if'i + 777 + azr]' + {2-5){l- 5)\^z-^ = n {2KZf] + r]^)"~^ [{2kz + 2r]){(p2 + \oz^~^) + 2K3Z^r]] . (40) 

In order to compensate the sub-leading singular behavior of /ij near x ^ we make the ansatz 

iP2iz) = -Xoz^~' + hiz) , (41) 

which we substitute into equation (^). For the function h we get the equation 

h" +^rj + azrj' = n{2Kzri + r]'^)"'^^ {{2KZ + 2rj)h + 2K3Z^T]) . (42) 

Since the function 77 is regular near z ~ 0, the solution h turns out to be regular near z = 0, too. Therefore, the function 
z 1-^ /i(l'Zl) is twice differentiable near a; = if we impose at z = the boundary condition 

h'{0) = . (43) 

At infinity we need that lim^^oo <P2(^) — 0- We therefore require that 

lim {h{z) - Xqz^-^) = . (44) 

z — >aD 

A proof of the following proposition can be found in the appendix. 

Proposition 6 For all n > 4, there exists a unique function h : R+ — > R that satisfies equation with the boundary 
conditions and (p^. Near z = 0, the function h is of the form 

h{z) ^ ho + h2Z^ + . . . , 

with certain coefficients Hq and ft,2, and for z large h is of the form 

/i(z) = Aqz^--^ + -— + .. . , 



with Xq as defined in (38), for a certain constant A', and with 5' as defined in Proposition 



4 The Reaction Front 

Using the properties of the functions fj,2, fi and (p2, we get from Theorem ^the following behavior of the reaction front 
F. 

Corollary 7 Let v be as in Theorem ^ and F as defined in (j7^. Then, for all z G R, 

|2 



lim t^-'^F{t'^z,t) = l{2n\z\rj{\z\) + rj{\z\rr = ^^"(Nl) = , 



V2 - 6??4 \z\ + ■■■ for |z| « , 

U2XKr/\zf + for \z\»l, 



and for all y =/= 0, 



i (2A.t)V + . • • for |y|«0, 



]imt-^F{V~ty,t)^li2^l, fi2n\y\)^{ 

I cxp(-n|y|V4)(2"-iC"/|2/r(^-'^) + ...) for \y\ » 1 

Here, rio, Va are as defined in PropositionUi and C is as defined in /Is^j. 



5 The Equation for tj) 



In order to simplify the notation we define the function u, 

u{x,t) = fii{\y\) , 

the function /^g, 



the function 



and the function 



M3(y) = ^2(2/) - ^2/ . 
{x,t) = t-''r^i\z\) + t-^^,M + t-'■^^,i\z\) 



4>i{x, t) = (j){x, t) — K 



Vt 



(45) 



(46) 



The function Wqo in Theorem |l| and in (|2^) can then be written as v^. = u + (j). 
Let now v — v^a + 'Z'- Then, 

{y^ - n'Y = ((u + + i,f - ti^)" = ((2u0 + <\?) + (2(u + 0) + ^) 

= E (2"<^ + '^')"~' (2(^ + 0) + V')' V'' • 

Therefore, substituting the ansatz v = + "0 into (|^) leads to the following equation for the function -0, 

i>^i\}" -V^-I- f (0) , 

with the function V , 

V = 2n {2u(j) + 0^)""^ (w + 0) , 

the function /, 

/= -0 + 0'/- (2U0 + 02)" , 

and the map T, 

f(V') =n (2^0 + 02) ""V + E (^) (2zI0 + 02)""'(2(u + 0)+iZ)'>'= . 

fe=2 ^ ^ 



(47) 
(48) 
(49) 

(50) 



5.1 The function V 

The function /i, //(y) = ^i{y)/y is strictly decreasing on R^., and therefore ^i{y)ly > Mi(^''y)/(^'^2/) for ^ > 1- Furthermore, 
the functions 77 and are strictly positive and is strictly increasing. These properties imply that, for i > t > 1, 
u{x,t) + t-'^r]{\z\) = ^i^{\y\) + t-^r{{\z\) > t-^(/Lti(|z|) + 77(|z|)) >t-'^co > 0, where cq = inf^>o(Aii(z) + ?y(z)). Next, since the 
functions fi^ and (/?2 are bounded and since 87 > e, for n > 4, we have that [^"^A'adyl) + t~^^(/'2(kl) | < const. and as a 
consequence (u + 0) and {2u + 0) are positive functions of x for all fixed t > tq, if tq large enough. 

Proposition 8 For n odd, n > 5, there exists tq > 1, such that for all t > Tq the function V is positive. 



Proof. The function {2u4 



'"^)" ^ is positive, for n odd. 



As a consequence, for n odd, equation (^) is of the form indicated in Section 2, with V — V and T ^ T. The rest of 
this section treats the case of n even, which, as indicated in the introduction, is slightly more delicate. It can be skipped in 
a first reading or if the reader is only interested in the case of n odd. 

So let n be even. The idea is to split V into its positive part V = V+ and its negative part Vi — V-^ , and to show that 
V- is small enough so that it can be treated together with the nonlinear term. Consider the function defined in (psf). The 
problem is that becomes negative for large values of x, and that therefore V becomes negative for large values of x. To 
understand why becomes negative, we note that the leading order term t^^ \z~^ in the large z asymptotics of t~'^ri{z) is 
compensated by the leading order term ~t^^\y^^ in the large y asymptotics of t^^ij.^{y). The leading order of at a: large 
is therefore given by the second order term in the large z asymptotics of 77 and the leading term in the large z asymptotics 



of ■ The first of these terms is proportional to t 



, and the second one is proportional to t z"^' 



-s' — 



t-'y^ 



y z 



The 



corresponding proportionality constants Aoo and A' can be computed for n > 5 and turn out to be negative. For n — A,5 these 
constants can not be obtained from asymptotic expansions, but numerical results show that they are in fact also negative in 

f.ViPQP mQPQ AA/p r\r\ nnf. tippH n nrnnf nf fViiQ mimprir'nl farf. Viprmiep ttip fnllmsrincr nrn-nnQif.inn le al«r» nrtrrf^ni: fnr Tin«if.i-^7P 1/ 



Proposition 9 For n even, n > 4, there exists ti > 1, such that the function Vi, satisfies for all t > ti the bound 

sup \Vi{x,t)\ < const. . (51) 

Proof. The idea is to write (f> as tlie sum of a function (pQ tliat is positive and a function (j)^ that absorbs the asymptotic 
behavior at infinity. Since n^iy) ~ XoV^^^ for y smaU, with Aq > 0, there exists yo > such that n^dy]) > 0, for all \y\ < j/o- 
Let c > 0, to be chosen below, and let 9 be the Heaviside step function, i.e., 9{x) = 1 for x > Q, and 9{x) = for a; < 0. 
Then, we define the function (j)^ by the equation 

cl,^ix,t) = -ct-''9i\y\~yo)\yf\z\-'' , 

and we set (f>Q = cj) — 4>^,. In order to prove that is positive, for c large enough, we write 4>q = 4>^q^ + 4''^\ where 

4'\x,t) = t-'{,j{\z\) X\z\-'9i\y\ yo)) + t-'''v,{\z\) + c - yo)t-'''\z\'-'' , 

and 

<|>^^\x,t)=t-^{^,,{\y\) + X\y\-'9{\y\-yo)) . 

f2) ('2') ('21 

(j>Q is positive for \y\ > yo, since in this case 0q {x,t) — H2i\y\) > ^^'^ 4>o positive for \y\ < yo by definition of yo- 
Next we consider 0^,^^ For \z\ < t^yo we have that (l3o\x,t) = t-'fri{\z\) + 1'^''' (p^Hzl) . But t-'<-q{z) + t-^'^Lp^iz) > for aU 
z > 0, and all t > t, if t is sufficiently large, since > 0, since is bounded, and since |</32(-^)l < vi^) for z large enough. 
Finally, using the asymptotic properties of 77 and we see that ^g^^ > for \z\ > Pj/g if c is chosen large enough. 
We now estimate the function Vi . From the definition of (j)^ we get that 

10^ (x, i)| < const. t-^(^'+i) , 

and therefore, since (/)g is positive, we have the lower bound 

for some constant c < 0, from which (j5^) follows. ■ 
5.2 The function / 

Theorem 10 Let I be as defined in ^4^), and let n > 4. Then, there exists a constant cj > 0, such that for all t > 1, 

/ dx \I{Vtx,t)\<ci t-^^^^ . (52) 
^ — 00 

The function / is even, and it is therefore sufficient to bound it for x > 0. The strategy of the proof will be to rewrite 
the function / as a sum of functions of the form t~'^G{y)F{z), with a > 0, and with G and F functions with appropriate 
asymptotic behavior at zero and infinity. Each of the terms in the sum can then be estimated with the help of Lemma |l^ 
below. In order to keep the notation as simple as possible, we suppress in what follows the arguments of functions whenever 
there is no risk of confusion. 

Proposition 11 For x > 0, the function —I is of the form 

n p 8 
p=2 q=0 1=2 



whe 



'n — 1 



A3 = n 



^4 = E ^pj + t-^T2y'-Pit-^''T3 + t-''^^+^'^Ti + t-^^n + t-^^Te)P , 

A5 = -t-^-^^3^P2 + azv>2) , 

A, = -t-^^-^-^-'nTr'lz^-'v (Ao - /.g/"') , 

A'j = t-2n7+27-e^2,n-i(2^^^g _ 2K\oy^-^) - t''^''''+'^''-'' n{2K\z^-^)''-^ {2^l^^l^ - 2kAo2/3-*) 
As = t-2"^+^nri"-i2r7((^i - ny) ~ n^y^) - t-^'^^+''n{2K\z^-')'^-^2Xz-\^l, - ny - K^y^) , 



where 



and where 



T2{y, z) = 2(^1 (y) - Ky)y-^z^T]{z) + 2ii^{y)^.^{y) , 
Tsiy, z) = {2^z + 2r^{z)) ^3(2) + 2^iMv'~^ , 
T4(y,z) = 2(^1 (y) - ^y)y2-*-'^*'-2^2(^) , 

T5(y) = M3(y)' , 



Proof. In terms of the functions (|54|)-(|59|) we get that, for x > Q, 

2u(f> + 4,^ ^ t-^'^Ti + t-'T2 + t"^^T-i + t-^'^^+^'^Ti + t-^'n + t-^^'n , 
and therefore {2u(f> + (jy^)'^ = Bi + J2t=2 where 

and where yl2, ^3 and A4 are as defined above. Since — / = 4> ^ 4>i + (2u(j) + (f^^ , it remains to be shown that 

n p 8 
p=2 q=0 i=5 

Using the differential equations for /i2, 77 and ip2, we find that 

Bi - <!>'; + = i?i + i?2 + 5'3 + ^5 , 

where as defined above, where 

Ri = -t-"^ {-{2nXy'-'r + (2/ii(A2/-' + A^a))") , 

and where 

+ n(t-2TTi)"-i(i-'*T3) 
The functions i?i and R2 can be further decomposed as foUows 



(2kA)"(— + ^)z-^ + n{T^Y-^ {(^2 + Aoz'-*)(2kz + 27?) + 2K3;2'r,) 



PA 



n p 



■2 J 



where i?^ and i?2 are as defined above, and where 



51 ^ -nt-27(n-i)-e(2«Azi-'')"^i(2(Mi - Ky)\y-^ + 2mi^3) 

^2 = ni^27(n-l)-eyn-lj.^ ^ 



It remains to be shown that 

3 8 

2 — 1 -1—6 

but this follows using the definitions. ■ 



5.2.1 Proof of Theorem ^ 

In order to characterize the behavior of a function near zero and infinity we introduce the following family of vector spaces. 

Definition 12 Let p and q be two real numbers withp + q > 0. Then, we define V{p,q) to be the vector space of continuous 
functions F from R+ to R, for which the norm 

\\F\\^,^=sup\F{x)\{x-^ + x'>) 

is finite. 

Note that, if a function is in V{p, q), then it is also in V(p', q') for any pair of numbers p' < p, q' < q for which p' + q' > 0. 
Furthermore, if Fi is in V{pi,qi), and F2 is in V(pi, 92), then the product F1F2 is in V{pi + P2,qi + 92)- 
The following Lemma provides the tool that we use to estimate the terms on the right hand side of (p3|). 

Lemma 13 Let F e V(^o: ^1) o.'i^d G E V(Go, Gi), and assume that 

1 - Gi < Fi , (60) 
1 + Go > -^^0 , (61) 

and that 

Fi^l + Go . (62) 
Then, there is a constant G > 0, such that for all t > 1, 

\G{x)F{t''x)\ dx < Ct-i , (63) 





where 

e-7-niin{Fi,l + Go} . (64) 

Proof. From ( |62|) it follows that either Fi < 1 + Gq or Fi > 1 + Gq. In the first case we get using ( |60|) that 1 — Gi < 
Fi < 1 + Go , and therefore 

/ \G{x)F{fx)\ dx< (snpx^' iFif^x)]) f -^\G{x)\ dx 

Jo \x>0 J Jo X 

< const, t^'^-^i , 

and in the second case we get using ( |6l[ ) that — Fo < 1 + Go < ^i, and therefore 

/ \G(x)F(t''x)\ dx< ( sup4^ |G(a;)| ) / x^°\F{t''x)\ dx 

Jo \x>0 J Jo 

< const. t-T(i+Go) _ 

■ 

We now show that the right hand side in ( |5^ ) can be bounded by a sum of terms of the form t~'^ G{y)F{z). For each of 
these terms we then show that the corresponding functions G and F satisfy the hypothesis of Lemma This then implies 
that t~'^ G{x)F{t'^ x) dx < const. t~^'^^^\ for a certain ^ depending on F and G. It is therefore sufhcient to prove that 
(T + f > 1 + 47 for all these terms in order to prove the inequality (|5^). 

Proposition 14 For y, z > we have the bounds 

T2{y,z) <f2{y) = 2\fi,- Ky\y-' (snp\z'v{z)\] +2\fi,fi,\ , 

\z>0 J 

Tsiy, z) < %{z) = \{2kz + 2r7) + 2 fsup \^i:i{y)y'^''\] z^^'tj , 

\y>o J 

n{y, z) < nAy) ■ Ti^z) = (2 \f,,{y) - Ky\ y^-'') ■ (z^'-^ \^,{z)\^ , (65) 

(5-21 \ 2-S I 



r6(y,z) <r6(z) = 2 sup A*3(y)y H 1^21 + 1^2! , 

\y>o J 

and Ti G V(0, (5 - 1), G V(3 - (5, (5 - 1), G V(2 - 5, (5' - 3), T4,i G V(5 - 5', <5' - 3), r4^2 e V((5' - 5, 0), T5 G V(4 - 25, 2(5), 
T^G V(4-2(5,(5' + (5-4). 

Proof. The inequalities (pSj) follow by using the triangle inequality and the asymptotic properties of the functions [i-^, 



Bound on the function A2. We have the bound 

n-l 

1^2 1 < const. ^ t-'^f^'' (rr'-Pfa 



p=i 

where a = 1 + 7 + 3pe^. The function G = T2' is in V((3 — Sep), and the function F = T"^^^^T3 is in V(2 — 5, 
3e(n — 1 — p) — 3 + 5'). Since 5' > 1, the inequahtics ( |60| ) and ( |6ll ) are satisfied and, since 5' < 5 for n > 3, 

r y-3£p ifp>2 ^„ „ 
Therefore, cr + ^> I + 7 + 3pe7 + (3 — 3pe)7 = 1 + 47 as required. 

Bound on the function A3. We have that = t^'^^^+'^'^Sa^+t-^^Sa^s+t-^TBg^g, where Bg^, = n{t-'^^Ti+t-^T2T~'^T,, 
z = 4, . . . , 6. Since \T2/Ti\ < const, t^^'^ , and e — 376 ~ 2-f, we have the bound 

1^3,41 < const, (Ti"-lr4,2) T4,l , 

with cr = 27(n - 1) + 7(1 + 5'). The function G = T^^i is in V(5 - 5', 5' - 3) and the function F = r"^^T4,2 is in V{S' -S,3). 
Since S' > 3 the inequahties (60) and ( |6l| ) are satisfied and 1^/7 = 6 — (5'. Therefore cr + ^= 1 — 37 + 7(1 + 5') + (6 — S')j = 
1 + 47 as required. 

Similarly, we have that 

t-^'\B3r^\ < const. t-^Tl'-' T5 , 

with cr = 27(n - 1) + 2e. The function G = T5 is in V(4 - 2S, 26) and the function F = T^""^ is in V(0, 3). The inequahties 
(|60| ) and ( |6l|) are satisfied and ^/-f — 5 — 26 — 3 — 6e. Therefore, a + £, ~ 1 — 37 + 2£+(3 — 6e)j — 1 + 47 as required. 
Finally, 

i"'^ 1^3,6 1 < const. t-''T^~% , 

where cr = 2(71-1)7 + 67. The function G = 1 is in V(0,0), and the function F = Tl'~^% is in V(4-2^, 3(n- l)£ + 5' + (5-4). 
The inequalities (|60| ) and (61) are satisfied and C/7 = 1. Therefore, a + ^ = 1 — 37 + 67 + 7= l + 47as required. 

Bound on the function yl4. Since the functions T3/T1 and Te/Ti are bounded, T4/T1 < const, t^^f and T5/T1 < 
const, t^'^'^ we have that 

n 

IA4I < const. ^t-27(n+p)yn < ^^^^^ ^~aj.n ^ 

where cr = 2717 + 47. The function G = 1 is in V(0, 0), and the function F — T" is in V(0, 3ne). The inequalities ( |6C| ) and 
(|6l|) are satisfied and ^/7 = 1. Therefore, cr + ^ = (1 — 7) + 47 + 7 = 1 + 47 as required. 



Bound on the function Ac,. We have the bound 

lAI <i""|37<^2 + «V2l , (66) 

where a = 2717 + 47. The function G = 1 is in V(0, 0), and the function F — \3^ip2 + cizLp'r^ is in V(2 — 6,6' — 2). The 
inequalities ( |60[ ) and ( |6l| ) are satisfied and ^/7 = 1. Therefore, 0- + ^= (1— 7) +47 + 7= l + 47as required. 

Bound on the function Ag- We have the bound 

\A^\ < const, t-'' (Tr^z^-'rj) |Ao - May'"'! , 

where cr = 1 + 7. The function G = |Ao - Msy'^^^l is in V(2, 2) and the function F = T{''^z^-^'q is in V(2 - (5, 1 + 2^). The 
inequalities (pO|) and (pi]) are satisfied and ^/7 = 3. Therefore, cr + ^= 1 + 7 + 37= l + 47as required. 



Bound on the function A-;. We have the bound 

\A-j\ < const. t~'"|Ti"-i - (2kAz1-'')"-1| |2/ii/i3 - 2K\oy^-^\ , 

where cr = 2717-27+5. The function G = |2;Ui^3-2kAoj/^"'^| is mV{5-6,6-3) and the function F = \Tl'^'^ - {2K\z'^-^y'-'^\ 
is in V(— 3, 3+(5' — 5). The inequalities ( |60| ) and ( |6l| ) are satisfied and ^/7 = 6 — 5. Therefore, cr+^ = (1 — 7) — 27+£+(6— (5)7 = 

1 4- ae rAmiirprl 



Bound on the function Ag. We have the bound 

l^sl < const. t-'"|ri"~i?7 - {2KX7^-^y'-'^Xz-^\ {fi^ - ny - Kgy^j , 

where cr = 2^7 - 7. The function G = - - ^sy^j is in V(5, -3) and the function F = \Tl'^^ii - {2k\z'^-^Y-'^\z-^\ is 
in V(— 3 — (5, 3 + 5'). The inequahties ( |6C| ) and ( |6l| ) are satisfied and = 6. Therefore, cr + ^= (1 — 7)— 7 + 67= I + 47 
as required. 

Bound on the functions A^ q. We have the bound 

< const, t-'^ KzSf"'' - \^~\2K\z^~'r~^\ (|2(/.i - Acy)^-^"' |2Mi/^3I') , 

where cr = 2n7 — 2p7 + pe. The function G = |2(/i]^ — Ky)y~^^ ^ \2iiiii^\'^ is in V(p(3 — 5), Sep + q) and the function 
F = |(z*77)P-'? r""P - A^'"''(2kz1-'^)"-p| is in V(-3e(n - p), 2 + ^' - 3pe). The inequahties @ and (H) are satisfied, and 

, f 5 - 3p£ ifp = 2 , 

' \2 + 5' -ipe ifp> 3 . 

Therefore, cr + C = 1 + 47, for p = 2 and (7 + ^=1 + 7(1 + 5') > I + 47, for p > 3, as required. 

This completes the proof of Theorem ^ ■ 

5.3 The Map T 

Equation (^) is of the form (j2^) if we define the map T by the equation 

I J (V'j + Vi'4) for n even , 

with T as defined in ( ^o|) and V^i as defined in Section |3.l| . Using the definitions, we see that T can be written as, 

nV') = EE^P.«V'''+', (68) 

p=l q=0 



with 



for (p, = (1, 0) and n odd, 

Vp^q = { Vi for (p, g) = (1,0) and n even, (69) 

Q)(^)(2?I<?!> + 02)«-p(2u + 2?i)P-9 forp + 9>2. 



Proposition 15 Lei Vp,q as in mdi). Then, for all i > 1, 



sup \Vp,q{x,t)\ < const. r^P'i') , (70) 



7(n-l)(<5' + l) /or (p,g) = (1,0) , 
e(p,g) = <( 27(n-2) + 27 /or (p, g) = (2, 0) , (71) 
27(ri-p) for {p,q) ^ (2,0) andp + q> 2 . 

Proof. The case (p, g) = (1, 0) follows from (pT|). Let now (p, g) 7^ (1, 0). Since e — 7 — 1^7 > 0, for all z^, < < (5, we 
find that 

supl^r |</>(x,t)| <i-^(sup Iz'' rKz)|+t-(^-^-''^) sup p-,iy)\ + r''' sup jz'^ ^2(^)1) 
xsR zeR+ aeR+ 2eR+ 

< const, i""*" . 
Furthermore, since fJ,i{y) — 0{y) near y = 0, 

1 2?I(^ + <i-^ 2^^4¥^ |z0| + 101' < const, t-^-^ . 



Since the function |u + 0| is bounded, it follows that \Vp^q{x,t)\ < const. t-<'P^'i\ with e{p,q) = 27(71 -p). For (p, 9) = (2,0) 
we improve this bound using additional properties of the function u + (p. Namely, since 2/(n — 2) < S — 1, we have that 



\V2,o{x,t)\ < const, sup |(2u0 + 

xeR 



< const, sup 

xGR 

< const, sup 



V y 



n-2 














V 



const, sup 



n-2 



< const. t-2("-2)7-27 . 



2r 



6 Proof of the main result 

For functions / in = Li(R)nLoo(R) we use the norms ||/||^ = / |/(a;)| dx, ||/||^ = sup^^^ |/(x)| and ||/|| = \\f\\^ + \\f\\^, 
and we denote by B the Banach space of functions ip in ioo([l, 00)) x for which the norm || 



m\B = ™P 



^47 



is finite. Let tq as in Proposition ^ and ti as in Proposition]^, and consider, for fixed r > maxjro, ti}, functions if) of the 
form 

^{x,t)=T-^^^{x/^,t/T) , 

with Lp E B. Let if be the fundamental solution of the differential operator dt — — tV{^x, rt), and let, for given v ^ J ^ 
the map be defined by the equation 



TZ{(p){x,t) = (fg^j^{x,t) + (fQ^2ix,t) +J\f{Lp){x,t) , 



where 



'fo.i{^^t)= K{x,t;y,l) iy{y) dy , 



R 



'^'0,2(2;,*) = T'*''r / ds / dy K{x,t;y,s) I{y/Ty,Ts) 



R 



and where 



AA((^)(a;,i) = ^^AAp,,(v5)(x,i) , 



with 



^p,q{^){x,t) = 



- ^47, 



/ dy K{x, t- y, s) Vp^.iV^y, ts) T-4^(f+«V(2/, 



R 



The integral equation ip = "^(((S) is equivalent to the differential equation (22) with initial condition ■0o(^) = ^{x,t) — 
t~'^'^v{x/^/t). We note that, since the function V is positive, the kernel K is bounded pointwise by the fundamental solution 
Kq of the heat equation. 



1 



1 



1 {x ~ yf 
4 {t-s) 



(72) 



Ko{x,t;y,s) = -7= ^- : exp 

The following Proposition makes Theorem ^ precise. 

Proposition 16 Let j3 > max|l, 3 cj ^1 + ^ ^ 1/2+4-7 1 ' wt/i c/ as defined in (|5^, anrf Zei r 6e sufficiently large. 
Then, for all v € J7 with ||z^|| < /3/6, t/ie equation ip = 'R-{ip) has a unique solution ip* in the ballh{{(3) = {ip €z B\ \\ip\\ig < /?}. 



Proof. Since 47 < 1/2, the solution ol the integral equation will be dominated by ipQ 2, and, as we wiU see, f3 has been 
chosen such that ||</'o2||g ^ /3/3. The idea is therefore to show that, if r is large enough to make the nonlinear part of 
the map TZ small, and if \\v\\ < /3/6, then the map TZ contracts the ball U{P) into itself, which by the contraction mapping 
principle implies the theorem. We first show that TZ maps the ball U{(3) into itself. For the contribution coming from the 
initial condition we have 



and therefore 



|(/.o,i|L <2|k|| </3/3 



We next estimate the norm ||<Po 2IL • Let c{t, s) ~ -j= 



Then, 



s/i \/t-s ' 

Vo,2{^ ■ ^t) <t"^^t dsc{t,s) I dy \l{y/Ty,Ts)\ 



R 



J \/lc{t,s)ds J dx \I{\/tsx^ts)\ 

< ci r^^rj' ^ c{t,s) ds {ts)-^^+*^^ 
ds 



3 



and therefore 



It remains to be shown that the nonlinearity is also bounded by /3/3, for r large enough. For ip £ U{P) we have, 



< 



c{t,s) \fs y~^( 

p=l q=0 



TS) '"'^^''^'^ s-47(p+9)^-47(P+9) ||^||^+'^ 



(73) 



For (p, g) = (1,0) we get, since 5i = 7(n - l){5' + 1) - 1 > 0, 

A/'i,o(¥')(V< ■ ,t) < const. T^^'l^ c(t,s) ds g-i-^'i-^T 

< const. r~*i/3 t''^'' 

and for (p, q) — (2, 0) we get 



' c{l,s) 

„l/2+47 



ds , 



AA2,o(<^)(Vt . ,t)j < const. ^47+l-87-(27(«-2)+27)^2y ^^^^ ^ ^_8^_2^_2^(„_2) 

1 c(l,s) 



,1/2+47 



< const. T-^P"^ I ?r„',7 ds 

and for the other cases we have 

Afp.q{ip)iVt . ,t) < const, ^47+l-27(r^-p)-47p-479^2„ / ^^^^ ^ ^-2^(n-p)-4jp-4jq 



< const. T-^^^" t 



7a2n ,-47 



' c(l,s) 

sl/2+47 



and therefore ||A/'(</3)||g < /3/3 if r is large enough. Using the triangle inequality we get that |l7?.((/?)||g < P, which proves 
that TZ {U{(3)) C U{[3) as claimed. We now show that TZ is Lipschitz. Let and (^2 be in U{l3). We have 

A/'(^i)(v^ . ,t) - M{if2){\^t . ,i)|| < const. ds c{t,s) ^fs ■ 

n p 

S-47(p+9)^-47(p+9) ||<^^-^2||g , 

and therefore we get, using the same estimates as for (|73|), that 



\\TZ{cp,) - 7^((^2)ll6 - - ^fi^2)\\B < 2 11^1 - ^2ll6 , 



iHprl T i« InrcTA pnmicrli TViiQ rnTnnlpf.pc; tTip r^rnnf nf Tlipnrpm 



7 Appendix 



7.1 Proof of Proposition ^ 

We first prove the existence of a unique positive solution of equation ( ^0| ) satisfying the boundary conditions ( |3T| ) and (|3^). 
Then, we derive the results on the asymptotic behavior near zero and infinity. 

7.1.1 Existence of the function 77 

Proposition 17 Let, for p > 0, rj^ be the solution of the initial value problem on H+, 

Tj" ^ {2KZf] + f]^y\ (74) 

r?(0) = p > . 

Then, there exists a unique p such that the function rj-^ is positive and satisfies linijj^oo Vpi^) — 0- 

Proof. We first prove that p is unique. Given a function r] from R+ to R we define the function T{ri), ^ iv) {^) — 
{kzt] + 77^)". Assume that there are two values Pi > P2 > 0; such that the functions rji = r/p^ and 772 = rj^^ are positive 
and satisfy lima;^oo ?yi(2;) = lima;^oo ^72(2;) = 0. We first show that the function r]i2 = Vi ~ V2 positive for all a; > 0. 
Namely, if we assume the contrary, then because 7?i2(0) > 0, there must be a first xq > such that 7?i2(a;o) = 0. Furthermore, 
if ryi2(a;) > then 7712(2;) = ^ iVi) i^) ~ ^ {V2) i^) > 0. and therefore 7712(2^0) = Pi - P2 + Iq" dx J^dy rfl^iv) > 0, a 
contradiction. Therefore rii2, and as a consequence 7712 , are positive for all x, from which it follows that limj;^oo Viii^) > 0, 
in contradiction with lima;^oo '?i2(^) = liin2;-+oo Vii^) ~ lim^^oo V2i^) — ^■ 

To prove the existence of a p for which 77-^ is positive and for which lim:r^oo Vpi^) ~ 0, we use the so called shooting 
method. Note that, for any p > 0, the initial value problem (|7j) has a unique solution r/^, and since 77p(0) = —k, the function 
77^ is strictly decreasing on [0, Xp) for Xp small enough. We will show that for small enough p > 0, the graph of ry^ intersects 
the real axis and 77^ becomes negative, whereas for p large enough, 77^ has a minimum and then diverges to plus infinity. The 
(unique) point between those two sets is p. Define the two subsets Ii and I2 of R+ , 

/i = {p G R+l 3 xi,f]p{xi) = and 77^(0;) > for x E [0, xi)} , 

I2 = {pe R+l 3 a;2,77;,(a;2) = and 77^(0;) < 0,77^(0;) > for x e [0,X2)} ■ 

We note that if 77^(0:0) = and 77^(^0) > 0, for some xq, then 77^ > on any interval (xq,x) on which 77^ is defined, and a 
function 77^ with p G I2 can therefore not converge to zero at infinity. Furthermore, since the function 77 = is a solution 
of the differential equation (|7^), it follows, since solutions are unique, that 77^(^0) > if ri'p{xo) — 0, and therefore the 
intersection of Ii with I2 is empty. The sets Ii and I2 are open, by continuity of the solution 77^ as a function of the initial 
data p. We now show that Ii is non empty and bounded, which shows that p = sup/i < 00. This p is neither in Ii nor in I2, 
and therefore the function 77^ is at the same time strictly positive and strictly decreasing, and therefore limj-^oo Vpi^) = 0- 
To prove that Ii is non empty, we fix any p^ positive and choose xq > small enough such that on [0, Xq] the solution 
r]i = rjp^ exists and is strictly decreasing. Then, p^ — 771 (xq) > 0. Choose now < P2 < Pi — ?7i(a;o) and let 772 = 77 be the 
corresponding solution. As before, we have that the function 77^2 = Vi^ V2j and its second derivative 7712, are positive on the 
interval [0, xo), and therefore, since 772(0:0) = P2 + /o^° !q '^V v'i.iy) = P2 + ^1(2^0) ~ Pi ~ /o"° dx dy r]'-[2{y), we find that 
772(^0) < P2 ^ Pi + '7i(2;o)- Using the definition of P2 we therefore find that 772(2:0) < 0. Therefore P2 £ h- We now prove 
that Ii is bounded. For p > 0, let Xp be the largest value (possibly infinite) such that on [0, Xp) the solution r/p exists and is 
strictly positive. Then, rjp = ^ {rjp) is positive on (0, Xp) and, therefore r/p(x) > p ~ kx ior x £ (0, Xp). As a consequence, 
if the function 77p exists on [0,p/k], then Xp > p/k. Using again that rip{x) > p — nx we then find that r)p{x) > p/2 for 
X G [0, p/Ik]^ and therefore T [rjp^ > (p/2)^" on [0, p/2k], which implies that 7/p(p/2) > — k + (p/2)^"+-'^, which is positive if 
p > 2k^/^"+^. Therefore ?7'(a;) must be equal to zero for some x < p/n. Any such p therefore belongs to l2- If the function 77p 
ceases to exist before x = p/n it must have been diverging to plus infinity for some x < p/k which again implies that r]'p{x) 
must have been equal to zero for some x < p/k, and the corresponding p is in 12- ■ 

7.1.2 Asymptotic behavior of the function rj 

The function 77 is regular at zero, and the coefficients of its Taylor series at zero can be computed recursively. We have 
77(0) = 77q > and 7;'(0) = ~k, and therefore we get using the differential equation that 7/2 = 7;"(0)/2 = 7;q"/2, ri"'{0) = 
and 774 = — 7;™(0)/4! — j^rjQ^^'^ {n"^ — tjq"^^). The asymptotic behavior of 7/ at infinity is obtained as follows. Assuming that 
77 behaves like X/ at infinity we get from the differential equation that S and A are as defined in (^6|) and (|3^), respectively. 
That this is indeed the correct leading behavior of 77 at infinity can now be proved by using standard techniques based on 
repeated applications of I'Hopital's rule. See for example [Q. Since the proof is simple, but lengthy and quite uninteresting, 
we do not give the details here. 

Once the leading behavior of rj at infinity has been established we make the ansatz ?7(z) = A^:^* + s{z). To leading order 
we get for the function s the linear equation 



s" - ji2K\rz-h = n{2>i\r-^\^z~^- 



(75) 



There is a certain constant Xp, such that the function Sp, Sp{z) = XpZ~^~^^ is a particular solution of equation rtTq). The 
solutions of the homogeneous equation associated with (pq) are of the form s^{z) — zP±, and using the definition ( |33D for A 
we find that 

P± = ^(l± Vl + 4n5((5 + l)) . (76) 

For n > 6 we have that \p-\ > 2S + I, and the asymptotic behavior of s is therefore for n > 6 of the form Xryo/z^^^^, with 
^oo = Ap, and of the form Aoo/^^'^"' with some unknown coefficient X^o for n < 5. It is tedious, but not difficult, to prove 
that this is indeed the correct second order behavior of rj at infinity. We omit the details. 



7.2 Proof of Proposition |5| 

In order to study equation (jsj) with boundary conditions ( |35|) and (|36|), we make the ansatz iJ.2{x) — m{x)/x^. For the 
function m we get the differential equation 

rn" + (^ - -W +(s{S + l)\-(^--e])m= l(2^m)" , (77) 



■2 X 

and the boundary conditions for m are 

lim m{x) = X , (78) 
lim m{x)x^^^^ = . (79) 



7.2.1 Asymptotic behavior of the function 

a solution of equation (F 



As indicated in Section p. 3 
exp(— a:^/4)/a;-^ 



The proof is similar to the one in 
we find for the function m at infinity either a behavior proportional to x^^^' 
Since 5 — 2e > 0, we find that 



I that is defined on R+ behaves at infinity either like x~^'^ or like 
We omit the details. Given the asymptotic behavior of /ij at infinity, 
or a behavior proportional to x^^ exp(— a;^/4). 



lim m{x) = 



(80) 



if and only if the boundary condition is satisfied, and we will impose ( pO| ) from now on. We now discuss the asymptotic 
behavior of the function m near zero. From equation ( |77| ) we see that m"(0) exists if and only if S{S + 1)to(0) = (Km(0))", 
i.e., if m(0) = A, and if to'(O) = 0. We then find, that m"(0)/2 — Aq, with Xq as defined in (38). By taking derivatives of 



equation ([7^) we find that to"'(0) — 0, and that to'^(0)/4! — Ai, for some constant Ai ^ 0. By taking further derivatives. 



one can recursively compute the Taylor coefficients of a solution toq of equation (|7^) that is regular (in fact, analytic) in a 
neighborhood of zero. The solution itiq does however not satisfy the boundary condition (^0[). The solution of ( [tt]) that does 
satisfy (^) is of the form 



m{x) = rriQ^x) + x^mi{x) 



(81) 



where p = p+ + with p+ as defined in (^6[). Here, mi(a;) = mi(0) + . . . , with m'i{Q) = 0, and with mi(0) to be determined. 
The asymptotic form (|8l| ) can be obtained by substituting the ansatz ( ^ ) for m into equation ([ttI). Since p > 7 we find 
from ( ^l| ) that near zero raQ{x) + x^mi(x) = A + Aoa;^ + Aix^ + . . . . We omit the details of the proof that the asymptotic 
behavior is as indicated. 



7.2.2 Existence of the function 

We now prove the existence of a function m that satisfies equation ( [TT] ) with the boundary conditions ( |78| ) and (pO[). Since 
the second derivative of the solution ni at zero is positive, and since m converges to zero at infinity, there must be a first 
^ G R^, such that m'(f ) — 0. The basic idea is now to use this position ^, and the value p of m at ^, as parameters in shooting 
arguments towards zero and infinity. The first shooting argument will allow us to define a curve cq of initial conditions (f , p) 
in R^, for which the boundary condition at zero is satisfied, and the second shooting argument will allow us to find on this 
curve an initial condition for which the boundary condition at infinity is satisfied as well. 

So, let (^,p) be an initial condition. Locally, i.e., near ^, there exists a solution m^^p of equation ([zTj). By definition, 
m^_p(^) = p, ^(^) = 0, and therefore we get for the second derivative of to^,p at ^, 

<p(e) = '^i(0p"+^2(c)p, 

where 




and 



■n 5(5 + 1) 

7r£ ^ 

2 f 



-^iO = ^e-"-^. (83) 



For initial conditions such that p — C2(^), where 



and — ^5{5 + 1)/ (^), we therefore have that »Ti^'^(^) = 0. See Fig. 1 for the graph of the function C2. The function C2 
has a maximum at the point that satisfies the equation 

<iUc2iW-'+io',{U = 0, (85) 

and the hue C2 divides the set of initial conditions into two subsets, a subset A where "^^'p(C) < 0, and a subset B where 
p(C) > 0. For initial conditions on C2 we can compute m^"^^^^.j (^), 



<'c.«)(0 -^(0C2(0"+^2(0C2(0 



Comparing with (^5|) we find that m^" ^^^-^ )('?m) ~ 0: ^-^^d we have that m^"^^(.^^(^) < for < ^ < i^„j. We now construct 
the fine Co for < ^ < i^„j. 

Proposition 18 Fix 5, < ^ < Then, there exists a unique number co(Oi C2(C) > co(C) > •^i sitc/i t/iat mj^co({) *'S 
positive and satisfies lim.x-^Q m^ co(i)i^) = Furthermore, the function cq is continuous. 

Proof. The proof is similar to the one in Section |7.1.1| . Define the two subsets Ii and I2 of the interval / = (A, C2{£,)), 

/i = {/9 e /| 3 < ^1 < m^,pi(^) = A and A < m^^p{x) < 02(2;) for x G (^i, 0} , 

/2 = {/9 e /| 3 < $2 < m^A^^) = 22(^2) and A < m.^Jx) < C2{x) for x G (^2,6} • 

The intersection of Ii with I2 is by definition empty, and the sets Ii and I2 are open, by continuity of the solution m^^p as a 
function of the initial data p. We now show that all p sufficiently close to A are in /i , and that all p sufficiently close to C2 {£,) 
are in I2. This implies that co(^) = sup/i < C2(^), and co(^) is neither in Ii nor in I2, and therefore the function 
satisfies A < mo{x) < m2{x) for all < a; < ^, and therefore lima;_^o ^tij.coCC) (2^) = since lima;_+o C2 (2;) = A. So let (^,p) be 
an initial condition. Then, m^^p satisfies the integral equation 

m^.pix) = P + J y p{z) ^{m^,p{z),z) dz , (86) 

where 

exp(zV4) 
Pi"^' ^ ^25 ' 

and where 



r2(s, z) — uj\ (2)3" +0^2(2)3 . 

r2(s,z) is strictly negative for < z < and s « A, and therefore we find, like in the proof in Section^ that any solution 
with an initial condition p sufficiently close to A will cross the line m = A. Similarly, for an initial condition (^,p) close to 
(^,C2(^)) we can use that i7(s,z) « 0, and that dS^{c2(z),z) is strictly negative to show that the corresponding solution 
will cross the line C2. This completes the proof of the existence of co(^). To prove uniqueness it is sufficient to use that 
ds^i^s, z) > for (s, z) in the set C (see Fig. 1), and to integrate the difference of two solutions from their respective initial 
condition to zero, which leads to a contradiction, since both solutions have to be equal to A at zero. Finally, that cq is a 
continuous function follows from the continuity of m^^p as a function of p and ^ using the uniqueness of co(^). ■ 

We now prove with a second shooting argument that solutions with initial conditions ($,co(^)), with ^ w 0, become 
negative somewhere in the interval (^, 2), and that solutions with initial conditions (^, co(^)), with ^ w stay positive and 
diverge to plus infinity. 

Proposition 19 There exists a unique initial condition (^*, Co(^*)) such that the corresponding solution is positive 

and satisfies limj.^00 f^^' ,co{oi^) — 0- 



Proof. Define the two subsets Ii and I2 of the interval / — (0, 



Ji = U e /| 3 $1 > m^^c„(^)(^i) = and m^^^^,(^^){x) > 0, m^^^^(^)(x) < for a; G (^,Ci)} , 
/2 = e /| 3 ^2 > m'^,c„{o(^2) = and m5,,„(5)(a;) > 0, < for a; G (^,^2)} • 

By definition, the intersection of Ii with I2 is empty, and the sets Ii and I2 are open, by continuity of the solution ^0(5) as 
a function of the initial data ^. We now show that all ^ sufficiently close to are in /i, and that all ^ sufficiently close to 
are in I2. This implies that ^* = sup/i < is neither in Ii nor in I2, and therefore the function m^* co(5*) is positive and 
decreasing for a; > 1^* which implies that \imx^oc'm^',co(^*)i^) — 0- i^^t (C,co(^)) be an initial condition with < ^ < xo, 
with a:o ^ 1- The proof that such an initial condition is in Ii is rather lengthy and we therefore do not give the details 
here, but on a heuristic level it is easy to understand why such a solution is in Ii. Namely, near zero the asymptotics of the 
solution coC?) is "^5,co(5) (2^) = ^ + Aoa:^ + • • • + (mi)^ ^^^.^^ {0)xP + . . . , where p — p+ + S, and where (mi)^ ^^^^^^ (0) is such 
that co(C)(^) ~ Neglecting higher order terms we find that co(C)(^) ~ 2Ao^ + {mi)^ ^^^^^ {0)p(.^~^ , and we conclude 
that (mi)''^_^^(^) (0) « - (2Ao/p)/e-^ Therefore, (^^^-'/'^Z^^-'^) « -2Xo/e^^ « and m^^^^^^^^i^^^''/^^^^''^) « A, 

if ^ is small enough. Therefore, since "^^'co(5)(^) < for all < a; ^ 1, we find that m^_co({)(^) < co(^) — 2Ao(a; — 

and therefore mj_co(C)(2^) = for some a: < co(OC^^'^/ (2Ao) + ^, as claimed. Next, let {£,, Co(C)) be an initial condition with 
£,„i — xi < ^ < Cmi with < a;i < 1 to be chosen below. We now show that such an initial condition is in J2. For all 
Co > > C we have the lower bounds m^^^g^^-j{£,') > 771^^^0(5) (C)^i and P(C')"i5,co(?) (^') - "^?,co(?) (0^2, where 



exp(r/4) Co " 



/ci > 0, and 



«0 



1 1 



k2 = (5((5 + 1) exp(-e^/4) / z-^' - ^) dz 



and for a; > we therefore have the lower bound 

m^^^„(^){x) > m5,co(?)(0 \ki+k2 



dy \ 

and it follows, using again the integral equation (^6|), that co(4) diverges at (or before) infinity, that therefore co(4)(^) ~ 
for some x > ^, which implies that ^ £ I2, provided 

fci + fc2 4^ > . (87) 

For xi small enough and for n large enough ( ^7|) can be verified without too much difficulty. With the help of a computer 
one can show that ( ^7|) is satisfied for the remaining n > 5. For n — i (^) is not satisfied, since the above bounds on 
™4 co(4)(Co) ^'^d co({)('Co) ^^^6 ^00 weak. Sufficiently good bounds can be obtained by dividing the interval (C^g) in two 
pieces and by integrating lower bounds on each of the subintervals. We omit the details. Finally, uniqueness of C* can be 
proved by integrating the difference of two solutions from to infinity, which, using the positivity of ds^{s, z), leads to a 
contradiction with the fact that both of the solutions converge to zero at infinity. ■ 

7.3 Proof of Proposition |B| 

We first proof the existence of a unique solution of equation (^2|) with the boundary conditions ( ^ ) and (^^. Then, we 
derive the results on the asymptotic behavior near zero and infinity. 

7.3.1 Existence of the function (p2 

The equation (^) for h is linear. We therefore first construct two linearly independent solutions hi and h2 for the corre- 
sponding homogeneous equation, which we then use to construct, using standard methods, a solution of ( p^ ) that satisfies 
the boundary conditions ( p3| ) and (|4^). The homogeneous equation is 

h" -qh^O , (88) 

where 

q{z) = n {2kz7j{z) + 77(2)2)"^^ {2kz + 2r]{z)) . (89) 
Since the equation (|88| ) is linear, the integral equation for hi, 

ry 

hi{x) = l+ dy q{z) hi{z) dz , (90) 



has a positive solution that exists for all x in R+. By definition, we have near x — Q the behavior hi{x) = 1 + O(x^). At 
infinity, the solution hi is asymptotic to a solution of the equation 



n 1 
/i"(x)--(2kA)"^/i(x) =0 



This equation is the same as the homogeneous part of equation (|75D, and the leading order behavior of hi at infinity is 
therefore either proportional to xP+ or to x^^ , with p± as defined in (|7^). Since hi is positive, we find using (pO[), that 
hi{x) > 1 for all x in R_|_, and therefore hi is near infinity of the form hi{x) = dixP+ + . . . , for some constant di > 0. A 
second solution of the homogeneous equation (^) is 

h2{x) = hi{x) ■ 

Near x = we have that h2{x) = a; + . . . , and near infinity we find that 

h2{x)=hi{x){d~d2X^-^^+ + ...) , (91) 

where d = l//ii(?/)^ dy, and d2 = {l/dif' / {2p+ — 1) . We note that hih'2 — h[h2 = 1. Therefore, the function hp, 

hp{x) = ci{x) hi{x) +C2(x) h2{x) , 

where 



ci{x) = h2{y) f{y)dy , 
Jo 

C2(x) = / hi{y) f{y)dy , 



and where 

f{x) — —jri{x) — axrj' [x) + n (2Kxrj{x) + 77(2;)^)" ^ 2K^x'^rj{x) 



satisfies equation ([421). Near zero, the function / is of the form f{x) — —^rjf^ + . . . , and therefore, using the behavior of hi 
and /i2 near zero, we find that ci is near zero of order O(x^), and C2 is near zero of order 0{x). The function hp is therefore of 
order 0{x^) near zero. At infinity, the function / is of the form f{x) = fooX~^ + . . . , where foo ~ —■yX + aXS + n{2KX)"~^2K3X, 
and therefore the function C2 is near infinity of the form C2{x) — di/co x''+"*'^~'^/(p+ + 1 — (5) + . . . , and ci is near infinity 
of the form ci{x) = —d C2{x) + h^o + <^i<^2/oo x^^p-'-^^ / {2 — p+ — (5) + . . . , for some constant hoc- Using these asymptotic 
behavior for ci, hi, C2, and /i2, we find for the fimction hp near infinity the behavior, 

hp{x) = i-d C2(X) + h^ + ^1^^2/00 2-p^-^ ^ ^ ^ ^ ^^^(^^ ^ C2{x)hi{x) {d - d2x'-^P+ + . . .) 

= hM^) ^li^did2xP^+'-'xP+x'-^P^ + -^lM^dix^-P+-'xP^ + ... 

p+ + 1 - S 2-p+- 5 



haohi{x) + ( 



„2-<5 



2p+ - 1 \p+ + 1-5 2-P+- 
= h^hi{x) + Xqx^-^ + . . . . (92) 

In the last equality we have used the definition (|3^) for Xq. The function h, 

h(z) = hp(z) - hoo hi{z) , 

solves the equation (^), satisfies the boundary condition (^), and since, as we show in the next section, the higher order 
terms in (^2|) converge to zero at infinity, it also satisfies the boundary condition (^). 

7.3.2 Asymptotic behavior of the function (/jj 

By construction, the leading behavior of h at infinity is h{z) = X^z"^^^ + . . . . We therefore make the ansatz h{z) — 
Xqz^^^ + fc(z), and to leading order, we get for the function k the linear equation 

k" - n5{5 + l)z-'^k = CkZ-^' , (93) 
for a certain constant Ck- The general solution of equation (ra) is 



k{z) — + const. zP- + const. 2^+ 



with a certain constant A' and with p_ as defined in ([7^). Since lim^^oo k{z)/z^^^ = but > 2 — S, the coefficient 
of the term proportional to zP+ must be zero. Therefore, since|p_| > 6' — 2 for aU n, the asymptotic behavior of k is always 
given by X' /z^ We omit the details of the proof that this is indeed the correct second order behavior of k at infinity. 
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Fig. 1. Initial conditions for m^^p for n — 6. 



